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INTRODUCTION 


This  supplement  describes  a revised  version  of  the  computer  program,  now 
designated  as  ARPY  5,  which  was  developed  to  analyze  stress -strain  data.*  The 
revisions  to  AMMRC  TR  76-12  are  as  follows. 

1.  Inclusion  of  a subroutine  to  arrange  the  input  data  in  increasing  order 
of  stress -strain  magnitude.  (In  the  original  version  the  program  was  aborted  if 
the  input  data  were  not  ordered.) 

2.  Inclusion  of  a data  extrapolation  routine  for  the  case  where  the  stress- 
strain  values  of  the  last  data  point  are  just  below  those  of  the  yield  values. 
(In  the  original  version  if  the  last  data  point  was  not  equal  to  or  larger  than 
Si,  the  Ramberg-Osgood  secant  yield  stress,  the  calculations  were  aborted.) 

3.  Inclusion  of  a new  method  for  tangent  modulus  calculations. 

4.  Inclusion  of  computations  for  the  correlation  coefficient. 

5.  Inclusion  in  the  output  of  the  stress-strain  parameters  expressed  in 
S - I metric  units. 

6.  Elimination  of  all  plastic  strain  component  analyses. 

7.  Miscellaneous  revisions:  Elimination  of  optional  Table  2 (tangent  modu- 

lus) and  Table  3 (plastic  strain  parameters).  Inclusion  of  S-I  metric  unit  values 
in  optional  Table  1 (comparison  of  input  data  with  computed  stress  values) . Compu 
tations  of  the  elastic-plastic  values  of  Poisson!s  ratio  have  also  been  eliminated 

Appropriate  changes  in  the  output  formats  have  been  made  to  accommodate  the 
various  revisions.  There  have  also  been  some  small  changes  made  in  notation,  in- 
ternal computations,  input,  or  output.  A few  revisions  have  been  made  in  the 
input  data  requirements  but  these  are  limited  to  data  required  on  the  header  card. 

Descriptions  of  the  revisions  are  given  in  this  report.  Samples  of  the 
printed  output  and  listings  of  the  revised  main  program  and  the  subroutines  are 
given  in  appendixes . 


DATA  ORDERING 

The  routines  for  obtaining  the  Ramberg-Osgood  S1  and  m parameters  require 
that  the  data  be  in  increasing  order  of  stress-strain  magnitude.  The  program 
contains  a test  for  data  order.  If  one  or  more  data  cards  are  not  in  the  required 
order,  a subroutine,  called  ORDER,  orders  the  data.  The  subroutine  contains  a 
provision  for  aborting  the  ordering  process  if  the  number  of  iterations  exceeds 
the  number  of  data  points.  When  this  occurs  the  calculations  for  the  given  data 
set  are  aborted  and  a diagnostic  message  appears  in  the  printout. 


*PAPIRNO,  R.  Computer  Analysis  of  Stress-Strain  Data:  Program  Description  and  User  Instructions.  Army  Materials  and  Mechanics 
Research  Center.  AMMRC  TR  76-12,  April  1976. 


1 


DATA  EXTRAPOLATION  FOR  RAMBERG-OSGOOD  ANALYSES 

In  the  original  version  of  the  program,  the  stress  value  of  the  last  input 
data  point  S(N)  must  equal  or  exceed  the  secant  yield  strength  Si,  or  the  Ramberg- 
Osgood  and  subsequent  analyses  are  aborted.  Let  a test  stress  S(T)  be  defined  by 

S (T)  = 0 . 7E/E  (N)  (1) 

where  E(N)  is  the  strain  value  associated  with  S(N) . Now  let  a test  increment 
be  defined  by 

AS  = S (N)  - S (T) . (2) 

Data  Test 

If  AS  is  negative  or  zero,  the  last  data  point  will  have  a stress  value 
equal  to  or  greater  than  S^.  If  AS  is  positive,  the  case  shown  in  Figure  1,  then 
S(N)  will  always  be  less  than  Sj.  In  the  revised  program  an  additional  test  is 
done  if  the  latter  situation  occurs.  The  ratio  S(T)/S(N)  is  calculated.  If  the 
value  of  the  ratio  is  0.95  or  greater,  an  additional  data  point  is  extrapolated; 
if  the  ratio  value  is  less  than  0.95,  the  R-0  calculations  are  aborted.  The  value 
0,95  was  chosen  arbitrarily  to  be  the  criterion  of  a last  data  point  which  was 


Extrapolation  Procedure 

A curve  of  the  strain-hardening  type, 

S = AeP,  (3) 

is  fitted  to  the  last  four  data  points  of  the  set 
using  a logarithmic  least-squares  procedure.  The 
stress-strain  coordinates  of  the  intersection  of  this 
curve  and  the  secant  modulus  line,  0.7E,  are  now  cal- 
culated. This  point  is  designated  S^,  e^  in  Figure  2. 
The  strain  difference  between  e^  and  E(N),  designated 
Ae  in  Figure  2,  is  determined. 

Now  2Ae  is  added  to  E(N)  for  the  extrapolated 
value  of  strain.  The  extrapolated  value  of  stress  is 
then  calculated  from  Eq.  3,  using  the  extrapolated 
value  of  strain.  This  procedure  assures  that  the 
extrapolated  data  point  will  always  be  in  excess  of 
the  secant  yield  stress. 


The  program  now  treats  the  extrapolated  point 
as  if  it  were  the  last  data  input  point  for  R-0 
computations.  However,  a diagnostic  message  appears 
in  the  output  printout  advising  the  user  that  the  ex- 
trapolation procedure  has  been  employed.  The  message 
also  contains  the  stress  and  strain  values  of  the 


close  to  the  Si  stress. 


Figure  1.  Identification  of  stress  and 
strain  values  used  in  the  program  for  a 
test  with  the  extrapolation  criterion. 
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Figure  2.  Extrapolation  procedure  outline. 


extrapolated  point.  No  entries  are  made  for  the  extrapolated  point  in  the  op- 
tional output  table;  a message  is  printed  to  this  effect  at  the  end  of  the  table. 


TANGENT  MODULUS 

If  the  Ramberg-Osgood  expression  is  differentiated,  the  resulting  tangent 
modulus  expression  is  given  by 

dS/de  « Et  - E/[l  + (3m/7)  (S/Si)”1-1]  . (4) 

The  corresponding  expression  for  the  strain  hardening  relation  is 

dS/de  = Et  = ApeP_1.  (5) 

Expressing  Eq.  5 in  terms  of  stress, 

Et  = Ap  CSi/A) q (S/Si)^  (6) 

where  q = (p-l)/p. 

When  the  value  of  E^  for  stresses  up  to  is  found  from  Eq.  4 and  for  stresses 
beyond  from  Eq.  6,  there  is  likely  to  be  a discontinuity  at  the  stress 
value.  This  discontinuity  can  have  two  forms  as  shown  in  graphs  in  Figure  3. 
Such  a discontinuity  leads  to  ambiguous  results  in  calculations  involving  the 
tangent  modulus  for  stress  values  at  Sj  (for  example,  in  plastic  buckling). 
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Figure  3.  Two  possible  types  of  tangent  modulus 
discontinuities  at  the  secant  yield  stress  (schematic). 


Note,  however,  that  it  is  possible  for  the  two  tangent  modulus  curves  to 
intersect,  thereby  eliminating  the  discontinuity.  The  stress  value  at  the  inter- 
section has  here  been  designated  as  the  equitan  stress  Sec^.  This  stress  value 
becomes  one  additional  parameter  to  describe  the  stress-strain  properties. 

The  equitan  stress  is  found  in  the  program  by  solving  Eq.  4 and  Eq.  6 simul- 
taneously, using  an  iterative  procedure.  This  is  done  by  a subroutine  called 
ISOTAN.  At  the  equitan  stress, 

E/[l  + (3m/7)(S/S1)m"1]  = Ap(S1/A)q(S/S1)q.  (7) 

Rearranging  Eq.  7 results  in 

CS/Si) q + (3m/ 7)  (S/S1)q+m~1  - E/[Ap(S1/A)q]  =0.  (8) 

The  S value  in  Eq.  8 is  the  desired  equitan  stress. 

It  is  possible  for  there  to  be  no  real  solution  for  Eq.  8 or  for  solutions 
which  are  not  useful  for  computation.  The  program  contains  tests  which  identify 
these  cases  and  diagnostic  messages  appear  in  the  printout.  The  nonsolution  and 
nonuseful  solutions  for  which  tests  are  done  are  as  follows. 

1.  No  solution  can  be  found  in  100  iterations. 
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2.  The  computations  diverge,  resulting  in  numbers  too  large  for  the  computer. 
(A  running  size-check  of  key  calculated  quantities  is  made  in  a subroutine  called 
SIZECK.) 

3.  The  equation  has  an  imaginary  solution. 

4.  The  solution  stress  value  is  larger  than  the  stress  of  the  last  data  point. 

5.  The  solution  stress  value  is  less  than  the  proportional  limit. 

Shown  in  Figure  4 are  the  results  of  a computation  for  a titanium  alloy.  In 
this  case  the  equitan  stress  is  less  than  the  secant  yield  stress.  The  curve- 
fitting parameters  for  the  stress-strain  data  are  as  follows: 

E = 15.2  x 106  psi  (104.8  GPa) 

Si  = 165.57  ksi  (1141  MPa) 
m = 25.07 


p = 0.1683 

A = 333.6  ksi  (2300  MPa) 

The  computed  value  of  the  equitan  stress  is  162.36  ksi  (1119  MPa). 


Figure  4.  Computed  tangent  modulus  values  for  a 
titanium  alloy  using  analytical  expressions  derived 
from  differentiating  the  Ramberg-Osgood  and  strain 
hardening  equations. 
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Referring  to  Figure  4,  note  that  there  is  a discontinuity  in  the  tangent  mod- 
ulus if  it  is  computed  from  Eq.  4 to  Si  and  from  Eq.  6 beyond  Si  at  the  values 
designated  B and  C in  the  figure.  The  discontinuity  disappears  when  Eq.  4 is  used 
for  stress  values  to  the  equitan  stress  (A  in  the  figure)  and  Eq.  6 thereafter. 

It  should  be  noted,  however,  that  the  stress-strain  values  associated  with  the 
tangent  modulus  would  be  computed  from  the  R-0  curve  to  Si  and  from  the  S-H  curve 
beyond  S1# 

Applying  the  procedure  results  in  a rough  approximation  of  the  tangent  modu- 
lus and  is  suggested  only  for  economy  of  storage  of  parameters.  If  a better 
approximation  is  required,  then  other  more  accurate  methods  should  be  used  which 
operate  directly  on  the  stress-strain  data.  For  the  latter  case  the  stress-strain 
data  may  also  have  to  be  stored. 


Metric  Units 


MISCELLANEOUS  REVISIONS 


The  output  parameters  are  given  in  both  Standard  English  units  and  S-I  metric 
units.  For  the  S-I  units,  the  modulus  of  elasticity  is  given  in  GPa  and  all  other 
stress  values  in  MPa. 

Header  Card 


An  input  Poisson's  ratio  is  no  longer  required.  The  following  changes  from 
the  information  given  in  Table  1 of  the  Supplement  to  AMMRC  TR  76-12  (page  4)  are 
required 


Columns  1-43 
Columns  44-50 
Columns  51-62 
Columns  63-65 


No  change 

Specimen  Area,  Format  F7.4 
Date  of  Test,  Format  2A6 

Output  Option,  Format  13;  enter  1 in  column  65 

for  optional  table.  Leave  blank  for  standard  output. 


Optional  Table 


Calculations  of  point-for-point  tangent  modulus  and  elastic-plastic  Poisson's 
ratio  have  been  eliminated  as  well  as  calculations  of  plastic  components  of  strain. 
These  data  were  printed. in  optional  Tables  2 and  3 of  the  original  program  and 
these  have  been  eliminated.  In  the  single  optional  table  in  the  revised  program, 
a point-by-point  comparison  is  given  of  the  actual  input  data  stress  values  with 
those  calculated  from  the  approximate  analytic  expressions.  The  values  are  ex- 
pressed in  both  Standard  English  units  and  S-I  metric  units. 

Additional  Parameters 


The  maximum  strain  of  the  data  is  printed  out.  Computations  of  the  correla- 
tion coefficient  are  made  for  three  separate  data  regions:  elastic  region;  from 

the  proportional  limit  to  the  secant  yield  stress;  the  strain  hardening  region. 
The  values  appear  in  the  printout  and  serve,  with  the  standard  error  of  estimate, 
as  an  indicator  of  the  goodness-of-fit  of  the  analytic  values  to  the  test  data. 
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APPENDIX  A.  SAMPLES  OF  OUTPUT 


1.  Sample  of  printed  output  of  stress-strain  parameters  and  statistical 
data  for  a specimen  of  ESR  4340  steel  tested  in  tension. 

2.  Sample  of  tabular  data  for  the  same  specimen  in  which  the  observed 
stress  values  are  compared  with  those  computed  from  the  Ramberg-Osgood  and  strain- 
hardening analytic  approximations.  The  columns  labeled  "Observed  Strain"  and 
"Observed  Stress"  are  the  input  test  data  values. 


ANALYTIC  APPROXIMATION  OF  STRESS-STRAIN  PROPERTIES 
RAM8ERG-0SG0C0  A NO  STRAIN  HARCENINC  PARAMETERS 
09TAINE3  FROM  ANALYSES  OF  EXPERIMENTAL  OATA 

R AM8ER G -0 SGC00  EOLATION  EPS=SIG/E+« 1/7  1 1 Sl/E H SIG/S1 I **M 


STRAIN  HAS  OENINC  LAW 


SIC-A  <EPS»**P 


SPEC  NO  ESR1C1  A 3 A 0 STEEL  TEMP=  20  C STRAIN  RATE=  1.0-05  PER  SEC  TEST  OATE  12-CC-74  MAX  ETRAIN=  .C193C 

*****  STANOARO  ENCLISH  UNITS  ***** 


HOOULUS 

S-INTCPT 

E-INTCPT 

SPL 

SI 

EXPNT 

.1-PCT  SY 

.2-PCT  SY 

COEFF-A 

SXPNT 

SEOT 

ME3A-PSI 

PSI 

MU-IN/IN 

KSI 

KSI 

M 

KSI 

KSI 

KSI 

p 

KSI 

29.52 

-94. 

3. 

162. 20 

174  . 

92 

103.90 

173.35 

174.52 

.2181*03 

.04623 

174.06 

I 

INTERNATIONAL  METRIC 

UNITS  ****< 

MOOULUS 

S-INTCPT 

E-INTCPT 

SPL 

SI 

EXPNT 

.1-PCT  SY 

.2-PCT  SY 

COEFF-A 

EXPNT 

SEQT 

SPA 

MPA 

MU-H/M 

MPA 

MPA 

M 

MPA 

MPA 

MPA 

P 

MPA 

203.51 

-.65 

3. 

1118.36 

120  6 . 

Cl 

1C  3 . 90 

1195.24 

12C3.24 

.1504+04 

.34623 

12CD.14 

STANOARO 

ERROR 

OF  ESTIMATE 

ANO  CORRE 

LATION  COEFFICIENT 

STANOARO 

ENCLISH  UNITS 

S- 

I METRIC  UNITS 

CORR.  COEFF. 

ELAST  SEE  R- 

0 SEE 

S-H  SEE 

rLAST  SEE 

R-C  SEE 

S-H  SEE 

ELASTIC  R-0 

S-H 

KSI 

KSI 

KSI 

MPA 

MPA 

MPA 

R R 

R 

755 

.259 

.395 

5.21 

1.73 

2.72 

• 

9999  .9929 

-98C2 

COMPARISON  OF  COMPUTED  STRESS  VALUES 
FCR  SPECIMEN  NO.  ESR1C1 


C9  SERVE 0 VALUES 
-CC-74 


0 8SERVE0  AO JUS  TEO 
STRAIN  STRAIN 


IN/IN 

OR  M/M 

.OOOOOC 

-.000033 

.0CC354 

.000351 

. 0C3740 

.000737 

.C01C66 

. CC1C63 

.001418 

.001415 

• CC 176  8 

. 0C1765 

•002122 

.032113 

• CO  2 59  3 

. CC2590 

.002343 

.002845 

.0C320B 

• CC32C5 

.003572 

.003569 

.003927 

. CC3924 

.004231 

.004278 

.004629 

.004626 

.005306 

.035003 

.005362 

.005359 

.005301 

.035733 

.006500 

.006497 

.007300 

.006997 

.007500 

.007497 

.008100 

.038C97 

•009C00 

• CC83  97 

.011200 

.011197 

•012660 

.012857 

.014300 

.014297 

.015000 

.014997 

.019300 

.019297 

STANOARO  ENCLISH 


8SERVE0 

CCMPUTEO 

STRESS 

STRESS 

KSI 

KSI 

.03 

- .03 

1C. 54 

1C. 35 

21. C3 

21.75 

31  .61 

31.37 

42.15 

41.76 

52.6  9 

52.09 

63.23 

62.54 

73.77 

76.44 

84.30 

83.97 

94.84 

94.59 

105.33 

105.34 

115.92 

115.82 

126.45 

126.26 

136.99 

136.54 

147-53 

147  .66 

158.07 

158.17 

168.61 

163.01 

173.00 

172.63 

173.50 

173.54 

174.00 

174.14 

174.73 

174.66 

175.10 

175.41 

176.90 

177.19 

177.95 

178.33 

179.00 

179.21 

1 80  .06 

172. SC 

181.11 

131.71 

WITH  EXPERIMENTALLY 
OATE  CF  TEST  12 


I METRIC  UNITS 


COMPUTES 

STRESS 

STRESS 

OEVIATION 

MPA 

MPA 

-.65 

.65 

71.39 

1.28 

149.95 

-4  .60 

216.29 

1.65 

287.92 

2.69 

359.15 

4 .13 

431.1 9 

4 .76 

527.04 

-18.42 

573. 94 

2.23 

652  .26 

1.73 

726.28 

.29 

798.52 

.72 

870.56 

1.35 

941.38 

3.13 

1013.11 

-.92 

1C9C  .55 

-.73 

1165.30 

-2.73 

1190.26 

2.53 

1196.52 

-.28 

12CC.64 

-.95 

1204.24 

.90 

1209.41 

-2.13 

1221.70 

-2. 0 2 

1229.53 

-2.61 

1235.58 

-1.42 

1238.32 

3.15 

1252.83 

-4.12 

UNITS 

STRESS 

DEVIATION 

KSI 

.094 
.185 
- .668 
.240 
. 39C 
.600 
.691 
-2.671 
.332 
.246 
.042 
.104 
.195 
.454 

- .134 

-.102 

- .403 
.367 

- .040 
-.138 

.13? 

-.310 

-.293 

-.379 

- .206 
.457 

- .593 


S- 

CeSFRVEO 

STRESS 

MPA 

.00 
72.67 
145.34 
217.94 
230 .61 
363.28 
435.96 
508.63 

581.23 

653.90 
726  .57 

799.24 

871 .91 
944.51 

1C  17 .18 
1089.85 

1162.52 
1152.79 

1196.24 
1199.69 
1205.13 
12C7.27 
1219.68 

1226.52 
1234.16 
1241.47 
1248.71 
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APPENDIX  B.  PROGRAM  LISTING 


The  program  is  written  in  FORTRAN  IV  and  has  been  run  on  a UNIVAC  1106  com- 
puter. The  listing  includes  the  main  program  and  five  subroutines.  The  number 
of  punch-cards  for  the  program,  including  comment  cards,  is  as  follows: 


Main  program: 
Subroutine  ORDER: 
Subroutine  LSTSQ: 
Subroutine  RJA75: 
Subroutine  ISOTAN: 
Subroutine  SIZECK: 


472  cards 
28  cards 
20  cards 
28  cards 
83  cards 
13  cards 


TOTAL:  644  cards 


SrCPtlS  MAIN 

COMMON  SI  IDO J .E(100 I .EE <1001  .X (100) .YI10C1. CS(100». 

ISM! 100  I *CSK! 100  I *0VM! 100 l»DV! 100  I 
DIMENSION  ARRAY! 1 ) 

E3UIVALENCEISI1 I » ARRAY 1 1 II 
SMC=6.  394757 
READ  7 0 » N1 
11=0 

1 I1=I1-*1 

IF! I1.GT.N1 ISO  TO  992 
DO  2 1=1*1000 

2 ARRAY! I )=C 

C READ  IN  OF  DATA  AND  CONVERSION  TO  PSI  AND  IN/IN  UNITS 

READ  71»SN0*XMTL  • YM  TL  •TEMd*EOOT*N*N3*P'>L*AREA*DATE3*DATE1!*N,7 
IF  1N3.EQ.2  IGO  TO  4 
IF! N3.E0.3 I GO  TO  6 

C RAW  DATA  IN  KSI  AND  MICRO-IN/IN  * N3=l  * 

DO  3 1=1, N 

READ  72.S1II.E1II 

S(  I ) =3(  I ) *1000. 

E ( 1 1 =E  (I  1/(10. **E.  I 

3 EE(H=E(II 
GO  TC  8 

C RAW  DATA  IN  POUNDS  LOAD  AND  TN/IN  * N3=2  * 

4 DO  5 1=1, N 

READ  7 3 , S ( I I , E 1 1 ) 

S (I l=S (I  I /AREA 

5 EEC  I l=E!  1 1 
30  TO  8 

C RAW  DATA  IN  KSI  AND  IN/IN  * N3=3  * 

E DO  7 1=1, N 

READ  74, Sill, ECU 
S ( 1 1 =S (I  1*1000. 

7 EE!  1 1 =E ! 1 1 

C CHECK  THAT  DATA  ARE  ORDERED  IN  INCREASING  STRESS/STRAIN  M ACNI TUD r 
3 NOUT=C 

NL=N-1 
DO  9 1=1 , NL 

IF  ( S 1 1 I • LE • S ( I* 1 1 ICO  TO  2 
IFCSI II.GT.S(I*1 } I NO  UT=NC JT  *1 
IF  (E  (I  I.GT.EI 1+ 1 1 | N CUT=  NO  UT* 1 
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CONTINUE 

IF  INCUT. EQ.C  IGO  TO  1C 
C ORDERING  OF  DATA 
MOPDER=L 

CALL  0RQER(N#3*Ef EE  #N ORDER ) 
IF (HORDER.Ea.OI  GO  TO  1C 
EN=E ( N I 
NEND=C 
GO  TO  51 

1C  S FL -PPL ♦ 1 L OG 

MX  -0 
SE  3T=0 
NRUN-0 

11  NRUN=NRUN*1 

C ELASTIC  DATA  ANALYSIS 
1=0 

12  1=1*1 

I Ft  St  I) .GT.SPL  IGO  TO  13 
IF  tl.EQ.NIGO  TO  13 
GO  TO  12 

13  JPL  =1 
J= JFL-1 

IF ( J.GT.2 ICO  TO  14 
EN  = Et  N ) 

C INSUF ICI E NT  ELASTIC  DATA 
NENQ=1 
GO  t0  51 

14  DO  15  I=ltJ 
X tl!=E  (II 

15  YfII=SIII 
L 1 = 1 
LJ=J 

CALL  LSTSO fXt YtSL # CFfLIt LJ  I 
EMOD=  SL 

IF (NRUN.E3.2 IGC  TO  1G 
SI  OF  = 0F 

E10F=-S10F/EMCD 
PE0F  = E1 OF* 1 1 0 • I 
PS0F=S10F 
EOFF  = E1 OF 
GO  TO  17 
1G  SO  OF  = CF 

E20F=-S20F/EM0D 
E3FF=E2CF 
PS0F=S10F*S2CF 
E120F=E10F+E2CF 
PE0F=E120F«t 10. ♦*£. I 
17  DO  13  1 = 1 » N 
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13  E ( I l=E  (I  i-corr 
C 7-0  ANALYSIS 
EN=E  (N  t 
SCT1=.7  *EM0D 
SCT2=.85*EM0D 
CHK1  = E(  N 1 * S C T1 
IF  (CHKl.GE.Sl Nl  IGO  TC  19 
CHK3-CHK1/S l N I 
IF  (CHK3.SE. .95  IGO  T0  AOO 
C LAST  DATA  POINT  STRESS  VALUE  LESS  THAN  SI 
NEND=2 
GO  TO  33 

C EXTRAPOLATION  ROUTINE 
**00  LI=N-3 
L J=N 
MX=MX*1 

DC  «I01  I=LI»LJ 
X(I1=AL0G(E(I)  J 
‘JOl  YCI»=ALCG (SCI l I 

CALL  LSTSQ(X.Y.SL.OF.LI.LJI 

XPN7=SL 

COEFF=EXP< OF J 

XTRPE=2.*(SCT1/C0EFFI**( 1 . / ( X PN T-l . II -E ( N I 

X7RPS=C0EFF*XTRPE**XFNT 

N=N+1 

E (Nl=XTRPE 
S<  NI-XTRPS 

19  I = J 

20  1=1*1 
DEL2=SCT2*E(I I 

IF < DEL2-SC I J J20. 22.21 

21  SL2  = fSCI  l-Stl-l  ||/mi  l-EI  1-1  II 

B?=Sm-SL2*E(I| 

S2=SCT2*B2/(SCT2-SL2 I 

E2=S2/SCT2 

JS2=I-1 

GO  TC  2 4 

22  S2=S(I( 

E2=E( II 
JS?=I 

23  1=1*1 

2«i  0EL1=SCT1»E(I  I 

IFC DZL1-SC 1 1 123. 2G.25 

25  SL1  = (S (I I -S 1 1 -1  I l/CEIII-Et 1-1 II 
B1=S( I J-SL1 *EI I ) 

S1=SC71*B1/(SCT1-SL1 f 

El =S1 /S  CT1 
JS1 =1-1 
GO  TO  27 

26  S1=S(II 
El  =E  I I ) 

JS1=I 

27  CHK2=JS1-JS2 

IF (CHK2.GE.2 IGO  TO  23 
N£ND=3 

C INSUFFICIENT  DATA  IN  KNEE  RECION 
GO  TO  33 
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23  JP-JS1+1 

XNUM=ALOG<  7.*(EM0D*E2-S2)/(3#*S1)  ) 

XDEN=ALOG  <S2/S1 1 
XM^XNUM/XDEN 

SPL=S1 *<7. ♦EMCD/l  Z.*Zl*l 10***G. ill *♦ ll./XK I 

IF< MRUN.EG.l ICO  TO  11 

IF  < MX  *GT  « 0 I ICHK-1GG**E(  N-  1 I/L1  + -5 

IF < MX«  GT  «G ) XTRPE=XTRPE*E12CF 

IF ( MX • GT • G I EN~E( N-l I 

EPL-SPL/EMOD 

S YrSl ♦ <7. ♦EM0D*.C02/C 3.* SI  I !**C 1 ./XM I 
SY1=S1*< 7«*EM0D*.0G1/ <3. ♦El (1 ./XM) 

DO  23  IrJPLtN 
23  Y ( I ) - AL  0 G(  S < I ) ) 

NP-N- JP*1 

IF<  NP« G E • 3 ) CO  TO  31 
DC  3C  I - JP  t N 
S<I)=C 
E < 1 1 =0 
EE< I )-0 
3C  Y(I|:0 

N=JS1 

C TNSUFFICENT  S -H  DATA 
NE  ND 

GO  TO  33 

C STRAIN  HARDENING  ANALYSIS 
C COMPLETE  DATA  SET  FOP  NEN3-5 

31  NENDrS 

DO  32  I-JPfN 

32  X (Il-ALCC  C E < X I I 
Li- Jr 

L J-N 

CALL  LSTS3(Xt Yf SLtCFtLlfLJ) 

P ”SL 

ArSl/( El  * *P ) 

33  D x 1 -0 
DT2-C 
DT3-0 
S3  A R1 -0 
SB  A P2  *0 
S3 A R3-0 
DRlrG 
DR2-C 

D R3r0 

CCR1 =0 

C 0R2 ~0 

C0R3  -0 

JROrJSl-J 

JSHrN-JSl 

DC  Zk  I=l» J 

CS( I ) -E  < I ) * E M 0 D 

OVII  US(I  l-CS(I  I 

DTI -DTI +DV ( I J * *2  . 

Zk  SBARlzSBARl+Stl  I/J 
DO  75  I=lfJ 
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35  0R1=0R1*ISBAR1-S( I II **2. 

DV1=S3RT( DTI /J ) 

C0R  = 1 .-DT1/DR1 
IF(C0R.LT.01  CO  TO  3G 
C0R1=S3RT (COR  I 

36  IF<  NEND.LT.M  I GO  TO  M2 
DO  37  I=JPL.JS1 

Z-C  S ( 1-1 J/Sl 
B E =EI  (I  I *E  MCD /S 1 
CALL  RJA75 ( Z * XM»  33 ) 

ccri  i=z*si 

CV( 1 1 =G  « IJ-C3CI J 
0T2  = DT2*DV  (I  1**2. 

37  SBAR2=S3AP2*  Sdl/J^O 
DO  38  I-JPLrJSl 

33  DR2tDR2*(jBAR2-5(I) |**2. 

D V2=SQRT (0T2/JR0I 
C0R-1.-DT2/DP7 
IF(COR.LT.OI  GO  TO  33 
CCR2-SQRT( COR  I 

39  IF (NCND.LT.5 I GO  TO  42 

DC  MO  I=JP.N 
CE  (I  I = A*E (I  I **P 
DV< 1 1=3 ( I J-CS( II 
DT3=DT3*DV(II**2. 

MO  SB AR3=SBAR3+S(I l/JSH 

DO  Ml  I=JP»N 

Ml  DR3=DR3+(SBAR3-S(II 1**2. 

DV3-33RT ( DT  3/JSH  1 
COR=2  .-DT3/DR3 
Ir( CCR.LT.O  I GO  to  43 
C 0 R3 -S  ORT ( CCR I 
GO  TC  M3 

CaUITAN  STRESS  COMPUTATION 

NO  R-C  DATA  OR  NO  S-H  DATA  * NC  E3UITA  N STRESS  COMPUTATION 
M2  JT  OUT  rl 

GO  T0  M 8 
M3  JT  CUT  rO 

SE3T-0. 

SE  GTM-O. 

XMTrXM-1 • 

XNT=(P-1. I/P 

CVJ=EMOO/(  A*P*  ( Sl/A)  »*XNT> 

KW  -0 
W=1  . 

CALL  ISCTAN (MtXMT  »XNT*CW»KWI 
IF( KW.LE.100  ) GO  TO  ME 
IF  CKW.LE. 20C I GO  TC  M5 
IF( KW.LE.300  I GO  TC  MM 

C NO  SOLUTION  CAN  BE  FOUND  FOR  ISCTAN  EQUATION 
JT0UT=2 
GO  TO  MS 

C MAGNITUDE  OF  VARIABLES  IN  IS''  T A N CONFUTATION  TOO  LARGE 
MM  JTCUT=3 

GO  TO  M3 

C IMAGINARY  SOLUTION  FOR  ISOTAN  EOLATION 
MS  JTOU  T-M 

GO  TO  MB 
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SCOT -W*  SI 

IF ISE3T. LE.SPL 1 GO  TO  47 
IF! SE3T.LE. SI ) GO  TO  48 
EEQTr<SEaT/AI**(l./PI 
IF ( EE9T • LT. EN ) GO  TO  43 
C EQUITAN  STRESS  LARGER  THAN  STRESS  OF  LAST  DATA  POINT 
JTOUT  = 5 
SE  QT -0 
GO  TO  43 

C E SUIT  AN  STRESS  LESS  THAN  Sf'L 
47  JT  OUT -6 
SE  QT -0 

43  IFIN2.EQ.0)  GO  TO  50 

IF  INEND.  LT.2  I GC  TO  51 
IF! NEN0.LT.4 ) NPTS-J 
IF (NEND.EQ.4 I NPTS-JS1 
IF! NEND.GT.4 ) NPTS-N 
DO  49  Irl.NPTS 
SI I)-S( I ) /1 000  • 

S M ( I I -S  1 1 USMC 
CSI  I ) -C SI  D/1000  . 

CSMII I-CS II  I *SMC 
DVI IUDVI I ) /1000  . 

49  D VMII  l=DV  II  USMC 

50  EMCD-EMCD/ 110. * *6 • I 
EMODM-EMOD*  SMC 
PSOFN-PSOF ♦SMC/1000. 

DVI -DVI /100 0 • 

0 V1M-DV1*SMC 

IF! NEND.LT.4  ) GC  TO  51 

XTRPS-XTRPS/1000. 

XTRPSM-XTRPS*SMC 

S1-S1/1G00. 

SI M-Sl * SMC 
SPL -SPL/1 COO. 

SPL  M- SPL ♦SMC 
SY-SY/100Q* 

S Y M- S Y ♦SMC 
SY1 -SY1/1000. 

SY1M-SY1 *SM  0 
0 V2-DV2/100U. 

Q V2  M-  D V 2 *SM  C 
IF (NEND.LT. 5 I GO  TO  51 
A - A/1 OD  0 • 

A M- A*5  MC 
QV3-DV3/10DP  . 

D V3M=D V3*SMC 

IF  (JTOUT. GT.O)  GO  TO  51 
SEQT-SEOT/IOCO. 

SE3TM-SEQT  * SMC 

51  PRINT  100 
PRINT  101 

PRINT  102  tSNOtXMTL  tYMTL#  TEMP . ED OT t D ATE 3 »D ATE 4 t EN 
IF! NEND.GT.l ) GO  TO  53 
IF  I NEND.E3. 1 I GC  TO  52 
PRINT  10  3 * SNO 
GO  TO  1 
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51 


54 

55 


5G 


57 


53 


EU 


93  9 
7 G 

71 

72 

"*■  7 
i -J 

74 

100 


101 

102 


^r/r  : c 4 

i r i N c.  \ - v_  : • •'.  i v;  C » C r'  L 
X.  «!;"N0.L3.r-  > cc  TO  57 
: r in  ' 1 r * 

ir  (WE  L'1*  i i CO  Tr  54 
PRINT  ICG 
0 0 "C  Nr 
PRINT  107 
PRINT  IDS 

PRINT  1D3.lM03.IH03m.P30F.  r'50FM«P20F»nVl  .DV1M.C0R1 
IF (N2.EQ.0 I GC  TO  1 
CC  TO  59 
PRINT  110 

IF(M5f.G'r  "j  ) PRINT  1 1 1 , 1 CH  K.  X \R  P3  . X T RPS  M . XT  RPE 
PRINT  111 

PRINT  11  3.:H03.PSOF.PROF.  I^L.Il  .XM.3Y1  .SY.3V1 .0V2.CCR1.  CCT-2 
PRINT  114 

PRINT  115,EM00M»PS0FM.PICr»3PLM»31M.XM.SYlM.3YM»DVlM»nv2M.CCRl. 
1C0P2 

IF(N?.E3.0»  GO  TO  1 
30  TO  59 

Icl  JTO'JT.ZO  .0  ) GO  TO  53 
PPINT  HE 

IF ( JTOUT.EQ. 2 I PRIN’  117 
I"( JTOUT.EQ. 3 ) PRINT  113 
IF ( JTOUT.E 0. 4 I FRINT  119 
IF( JTOUT.EQ .5  ) PRINT  120 
IF  (JTOUT.EQ. 6 ) F R I N ’’  121 
P R I N T 122 

PRINT  123.EM0D.P5CF.FE0r .3  PL  » 31  » XM  » S Y 1 » S Y « A . F . SE3T 
PRINT  124 

PRINT  125 . EMODM.rSCFM.PCOr .SFLM.SIM.XM.SYIM.SYM.AM.F.SEGTM 
PRTNT  1 2 o 

PRINT  127.DV1.DV2.DV3.DV1M.DV2H.DV3H.C0R1.CCR2.CCP3 
Ic(  N2.  E 3 .C  I GO  TC  1 
PRINT  128.SN0.DATE3.DATE4 
DC  GO  I-1.NPT3 

PRINT  1 2 9 . EE  ( I I »E(  I)  .5(1  » . C5  ( 1 1 .DV  (II  »SK(I  I.C3H,  (I>  »DVH(I  I 

IF(MX.GT.G)  PRINT  1 3T 

IF (MX  .GT.G  > 30  TO  1 

IF( NEHD.E3.4  ) PRINT  131 

IF (NEND.LT.4 ) PRINT  132 

CO  TO  1 

STOP 

FORMAL  13  ) 

rOPMAT (5AC  »2I2.F7.2»F7.4  .2AE.I3  I 


FORMATt  F7.2.F7.D  > 

FORMAT  (F7.DtF7.5l 
F3RHMI  F7.?tF7.5  1 

FORMAT (lHItTH2t 'ANALYTIC  A Pro OX jh A TI C N CF  STRESS-STRAIN  FRCFF  RTIE  3 
l’t/t  T(»5t  ■RAM3ERE-35G03  3 AND  STRAIN  HARDENING  PA  R 4 HE  TERC*  t /•  T4  0. 
2*0BTAINED  FROM  ANALYSES  ZT  EXPERIMENTAL  DATA'./I 
F3RMAT(T33t,RAM3E*’C“DSG303  EQ,JATIi:N,*EX«tEPS=SIC/IF*lj/7K3l/E1<ST^ 
1/31  |**^'»//tT3ot*STRAIN  HARDEN  INC  LAW' *5Xi*$I6=AlErS  t / I 

FORMATt  TIE.  1 O'*  EC  NO  * » A6  . 2X  t 2 AS  » 2X  t * T'  Mp3  1 . AG  • * C 1 1 TX * 1 0 T RA IN  RATE 
;=  * , AG  » * rnt  SE  C*  « 2*  t ' TEST  DATE  ’ t 2 AE  t?  X t * MA  X S TR  A I Nr  * . F 7 . S t / t 
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1 D 3 FORMAT! / * T2  5 » 3 (5  1 1 H * > » 6 X ) » // » T 32  . ******  NO  ANALYSIS  WAS  MADE  OF  TH 
IE  DATA  FOR  SPEC. NO.  *. A6  ,***•♦*’,// . T25 .* T HE  DATA  CARDS  COULD  NOT 
29E  ARRANGED  IN  INCREASING  ORDER  OF  STRESS  MAGNITUDE,*.//.  T25»9(5!l 
3H*  » *£  X I I 

10*1  FORMAT!/. T25. 315  IlHOtSX)  .//. T 27  ******  * INSUFFICIENT  DATA  IN  THE  E 
1LASTIC  REGION  WERE  INPUT  TC  THE  PROGRAM.  *****»./ ,T2 A .» A VALID  LEA 
2 ST  SQUARES  ANALYSIS  COULD  NOT  3E  HADE  AND  THE  CALCULATIONS  WERE  AB 
3CRTED. ’t//tT25. 8! 5( 1H*  ». EX  IJ 

105  FORMAT! T10. ***** *THERE  WCRE  INSUFFICIENT  DATA  FOR  A VALID  R-C  ANAL 
1YSIS  * I 

106  FORMAT! 1H+.T71 .’ TOC  FEW  DATA  POINTS  IN  THE  KNEE  OF  THE  S-E  CURVE.* 

1*****./! 

107  FORMAT! 1 H*  * T71 * ' STRESS  VALUE  OF  LAST  DATA  POINT  LESS  THAN  SI  STRES 

IE. ******  t / l 

1 D 3 FORMAT!  T A»  1 . •♦•*•*  ONLY  ELASTIC  DATA  ANALYZED  *•*•*»♦//, T20 .* ELA S TT 

1C  MODULUS  ’STRESS  INTERCEFT 'tTEtf  'STRAIN  INTERCEPT*  iTS2»*STND 

2 ERROR  OF  ES  T IM  A TE1  i HO  3 . » CDRR  EL' 1 t . T2Q  . * MED  A-  *>S  I*  » T ZZ . ’ GPA  ’ . T4  4 . ’ 
jrSI’.TSM,  »MPA,.T£4.’MlCRCETRAIN'rTS7»  »K5I*»T37»*MPA*.'103.’CCEFF*. 

4 / 1 

103  F0RKAT!T20.F7.2.T3r.F7.2.T42»F7.0fT52.F7.2.TC5.F7.0.T04,F?.3.T34.F 

17.2»T107»FG  .*1  ) 

111!  FORMAT (T14 .’*♦**♦  T HERE  WERE  INSUFFICIENT  INPUT  DATA  BEYOND  THE  SI 
1 STRESS  FOR  A VALID  STRAIN  HARDENING  A NA LYS IS .**♦**»,/ » 

111  FOPMAT  !T14  , •**♦♦♦  THE  LAST  INPUT  CATA  POINT  HAD  A STRAIN  VALUE  WHI 
1C  H W AS  * i 13.’  PERCENT  OF  El.  AN  ADDITIONAL  POINT  WAS ’ , / . TIM t *EXTR 
2APCLATED  SC  THAT  AN  R-0  ANALYSIS  COULD  BE  MADE.  THE  STRESS  AND  3TP 

3 AIN  VALUES  FOR  THE  EXTRAPOLATED  DATA  ARE  * . /.  T1  <1  . F 7 . 2 . ’ KSI  (’.F3.2 
*1.’  MPA  ) AND*  *F9.6t*  IN/IN  ! M/M  }.*****»,/ | 

112  FORMAT  !T5C» ******  STANDARD  ENCLISH  UNITS  ***♦*’.//. 

1T11. *MODULUS*.T21.*S-INTCPT*.T3C.*E-INTCPT*  « T 4 2 » *SPL*«  T51«*  SI*. T53 
2. ’EXPNT’.TEE . *. 1-FCT  S Y • . T 77 . ’ . 2-P CT  SY * . TB3 . ’ ELAST  SEE ’ . T SO . ’ R-0 
3 SEE’ .T1 DO. ’ELAST*. Til  3. *R-C* ./. 

4T11 . ’MEGA-PSI*  * T24 . *PSI* .T70. *MU-IN/IN* .14  2. *KSI*  »T5l*  »KSI *.TF1. *M 

5 *. TG9. »K5I* . T30. *KSI* .T31 * * KSI*  » T101 . * KSI ’.Till. *3* .Til  3. *R* . /) 

113  F0RMAT!Tll»F7.2»T22»F7.Cl»T3C!«F7.ri»T4L»F7.2»7"43,F7.2»T'"3,r-'.2»TC7»F 
17.2»T73.F7.2.T33»F7.3.T3S.F7.3.T1C3.F6.4.T11G.F6.4) 

11*1  FOPMAT  (1HL.T*1E.  ******  S-I  INTERNATIONAL  METRIC  UNITS  ♦****»,//, 

1 Til. ’MODULUS* .T21 .* S-INTC°T  * » T 30 . * E- INTCPT* . T 42 » *F PL*  » T51 * * SI  * . T 53 
2. *EXPNT’.T66. •. 1-PCT  S Y ’ . T 77 . * . 2-P CT  S Y * . T S 3 . * E LAST  SEE ’ . T 90 . * R-C 
3 SEE ** T1 C 3. * EL  A ST ’.Til  9. *R-C’./. 

<!Tl3»,CrA*iT2M.*MPA*fT31»,KU-M/M*.T*l?.*MFA*tT51.,MFA,.TCl.*M’.T63.’ 
S HP  A’  . T 3C  . *MPA**T91»  'MPA'tTlCl.*HPA»,Tllli'R'iTUa.»P*./» 

115  F ORM  AT  tTl  1 tF  7 . 2 ■ T 21 1 f 7*2  »T  3C1 1 F7 .0  i T3  E i F & *2  *T  4 B »FS  * 2 »T  5 3 » F 7 .2  » TG E » F 

13.2*T77»F3*?*T38*r7  .2iTSS.F7  .2  . TirB.FE.!!  iTl  IS  *FG  .«i  t 
HE  FORMAT  fTlLit1*****  TO  CONFUTE  THE  7ANCENT  MCOULUS  USE  THE  R-C  CURVE 

1 TO  SI  AND  ’’HE  3-H  CURVE  3EYONO  31.  AN  CCUITAN  STRESS  C DUL*1  * » /«  TIP 
2. ’NOT  BE  DETERMINED  BECAUSE’! 

117  FORMAT! 1 H* » T36. *N0  SOLUTION  COULD  3E  FOUND  FOR  THE  ISCTAN  EQUATION 

2 . * * * * * •,/ > 

113  FORMA"*!  1 H*  . T3G  . ’ THE  COMPUTATIONS  INVOLVED  NUMBER  SITES  TOC  LARGE  F 
10R  ’HE  COMPUTER.*****’./! 

119  FORMAT!  1H>.  T3S. ’THE  ISOTAN  EQUATION  HAD  AN  IMAGINARY  SOL U TI ON. * * * ♦ 

1*  * ./  I 

1 DC  FORMAT!  1 H* . T3S  .*  THE  VALUE  FOUND  WAS  LARGEP  THAN  THE  MAX-STRESS  DAT 

1A  POINT. *****  * . /I 

121  FORMAT!  I H* .* THE  VALUE  FOUND  WAS  LESS  THAN  THE  PROPORTIONAL  LIMIT.* 

1*****,/! 
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122  FORMAT!  T50*  * *****  STANOARO  ENGLISH  UMTS  •••••’,//. 

IT  11  . ’MODULUS*  .T21  . • S- INT  CP  T ’ . T30  . * E -INTCPT  ’ . T4  2 . * SFL  ’ . T51 . ’S 1 * » T5 9 
2 . * EXPNT  * » TG  6 » * ,1-PCT  S Y * . T7 7 . ’ .2- PC T S Y * . T 9 3 . ’ C OEFF -A ’ . TIOOt ’ EXP NT 
3’ .TIOS.’SEOT* ./. 

<1  Til . ’MEGA-PS  I’.T24.  ’PS  I’.  T30.’  MU- IN/ IN’  .T42.  • KSI’  .T  51.  ’ KSI’  . TS1.  ’ M 
5 * » TE  9 • 'KSI * * T80  » * KSI ’ » T9C » *KSI*»T102»*r’.T110.  *KSI  ’ » / I 

123  FORMAT!  Tl 1 . F7 .2 . T22 . F7 . 0 . T 30 . F7 .0 .T40 . F7 .2 . T 4 3.  F7 .2  . T 53  . F7 . 2 1 TG7.  F 
17.2»T78»F7.2»T87»E9.4fT2S»F7.5»Tir3fF7.2l 

12  4 FORMAT!  1HP.T4G.’  *****  3-1  INTERNATIONAL  METRIC  J NIT  3 

1T11  • ’MODULUS*  *T21  . ’3- INTCPT*  . T 30  . * F- 1 NTCPT  ’ . T4  2 . • SFL ’ . T51 . *S 1 ’ . T5 C 
2.  ’EXPNT’ .To  St ’ .1-PCT  S Y • . T7 7 . ’ .2- PC T S Y ’ . T 3 8 . ’ C OEFF -A ’ . Tl 00. ’EXPNT 
3*»T109»*SE9T*»/» 

4T13. *3PA».T24. *MPA’.T31.’  MU-  M/  M*  . T4  2 . ’ MPA  * . T 51  . ’ MPA  • • T S 1 . ’ M ’ . TG  3.  ’ 
5MFA ’ »T8C» ’MPA’. T90 . ’ MP A ’ *T 1C2 • * P * . Tl 1 0 . ’ MP A • • / I 
125  FORMAT!  T 11 . F7 .2 1 T21 . F7 . 2. T 3C . F7 .0  * T 39. F3 .2 . T4 3 . F 3 .2 .T S3 1 F7 . 2 . TGG . F 
iB.2»T77»F8.2»T87tE9.4»TS5»F7.5»T107.F8.ZI 

128  FORMAT!  1H0.  T40. ’STANDARD  ERROR  OF  ESTIMATE  AND  CORRELATION  COEFFIC 
1ICNT’.//,T20.  ’STANDARD  ENCLISH  UNITS  *.  T £ *» . *S-I  METRIC  UNITS*.  T99.  • 
2C0RR.  COEFF.’./t  TIS.’ELAST  SEE ’ • T2 7 . ’ R- 0 S EE ’ . T 39 . ’ S-H  SEE’.T5E.’ 
3ELAST  SEE*  .T63.  *P-0  SEE’  »T  8C»  ’S-H  SEE  ’ « T 94  • • E L AST  IC  ’ • T 1 C4  • ’ R-C  ’ . T 1 
412.  ’S-H*  ./.  T19.  ’KSI*  .T29.  ’KSI*  .T41 . ’ KSI*  .T53.  ’MPA’  .T70.  ’ MPA  ’ . T32.  ’ 
5MPA’*T97»*R’»T105»*R*»T113»*R’»/I 

127  FORMAT! T 1 5 • F7 . 3 . T27 . F7 . 3. T 33 . F 7 . 3 . T 5 7 . F 7 . 2 . TS 3 . F 7 . 2 . T 30 . F7 . ? . T34.F 
16  .4  »T102  » TG.  4 .T110.FG.4 I 

129  FORMAT 11H1.T27. ’COMPARISON  CF  COMPUTED  STRESS  VALUES  WITH  EX PF RIME 
1 NT  ALL  Y OBSERVED  VAl  UES’  »/.T3G»  * F O'5  SPECIMEN  NC.  ’ . A 5 . 3X  , • DA  TE  OF  T 
2EST  * «2A6  »//» 

3T47. ’STANDARD  EN  CL  IS  H UN  ITS  ’ . T 33  . • S- I METRIC  UNITS ’./.  T24  ,’  C3  S£R  VE 
40  *.T34 . ’ADJUSTED’ .T44. ’CBEERVED  *.T54.  *C OMF UTEO  * . TG5 • * 5 TRESS  * • T 7 S • ’ 
S03SERVED  ’ »T33»  ’ COMP’JT  ED’  . T92  * ’ STR  ESS  • . / . T25  . • ST  R A IN’  . ” 3 5 . • S TRA  IN  • . 
GT4  5 . ’STRESS’ »T5 5. ’STRESS ’ . TE4 . ’ DEV IAT I ON • . T73 . ’ ST RES S ’ . T 39 . * 3 T RE S S 
7 ’ . T33. ’ DEVIATION ’ ./. T25. • IN/ IN  OR  M/M’ » T47 . ’ KS I ’ »T 57 » ’ KSI* . TG7. »K 
ESI  * .T81. ’MPA’ .T91 . ’MPA’ . Tl 01. ’MPA’. /I 
123  FORMAT!  T24,F3.G.T34.F3.6.T44.r7.2.T54.'-7.2.TG4.F7.3.T73.F3.2.:3  3.F 
18.2  »T9S  »F7.2 I 

130  FORMAT! 1HC.T24.’ *****  NO  ENTRIES  FOR  THE  EXTRAPOLATED  DATA  POINT  A 
IRE  INCtUDED  IN  THE  TABLE  **♦**’! 

131  FORMAT! 1 HO . T24 .’ *****  NO  ENTRIES  FOR  DATA  3EY0ND  31  ARE  INCLUDED  I 
IN  THE  TABLE  *♦***’  I 

132  FORMAT! 1H0. T24  »*  *****  NO  ENTRIES  FOR  DATA  3EY0N3  THE  INPUT  VALUE  0 
IF  SPL  ARE  INCLUDED  IN  THE  TABLE  •****’! 

END 

BTORtlS  ORDER 

SUBROUTINE  ORDER <N. S. E. EE . MCRDER I 
DIMENSION  S! 1 I.E! 1 I.EEtl I 
COUNT-C 
MrN-1 

1 NOUT=C 

DO  4 I =1.M 

IF!  S!  II-S!I*1 114  .”>.3 

2 IF  !E  !I  I.LE.EII+1I  I GC  TC  4 
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3 NOUT - NOU  T+ 1 

S A -S  < I I 

EArE< I) 

S3=S<I + 1 | 

EG”E< 1+1 J 

S <1 IrSB 

E( IJ-E3 

S I I+l I -S A 

E( I+l I =E A 

CONTINUE 

COUNT-COUNT* 1 

IF  <COUNT.LT.NI  GO  TC  5 

MORDERrl 

RETURN 

5 I r ( NOUT  « GT  • C ) GO  TO  1 
DO  6 I -1  • N 

6 EE(  IJ=E( I) 

PE  TURN 
END 

2FC  ^ • IS  L STS  0 

SUBROUTINE  L ST  SI ( X t Y t SL . 0 F . L It  L J ) 
DIMENSION  X ( 1 I * Y( 1 ) 

ETrp 
STrp 
EET  -0 
E ST  -0 

DO  1 I-L It L J 
E T -E  T +X  (I  I 
ST-ST+Y  i I J 
EETrEET+X IX  I * ♦ Z 
1 EST-EST+X  1 1 1 * Y I 1 1 

C-LJ-LI+l 
DEN=C*EET-ET»*2 
ENUM-C*E3T-ET*3T 
CNUMrOT*EET-ET*EST 
Sl-ZN’JM/DEN 
OFrCNUM/DEN 
RETURN 
END 

SFCR.IS  RJA75 

SUBROUTINE  R JA7 51 2 tXMi B3  I 
FUNGI Z)-ZM3./7.!*Z**IXM1-BB 
KK-P 

RM=.OOOCI 

L=0 

3Z-Z/1GC  • 

1 RZrFUNCtZI 

I F ( A 3 G ( RZJ.LE.RM ) C C TO  7 
KK=KK+1 

Ir< KK.GT.lOniGO  TO  7 
Ir  IRZ.LT.U.  I G C TC  3 

2 L=L+1 
Z=Z-L*DZ 
GO  TO  1 


17 


3 V=Z* OZ 

RV=FUNC<  V) 

IF  ( ABS  ( R V I.LE.RMI CC  TC  6 
KK-KK+1 

IF  I KK-IOP  ? 4 • *1  • 7 

<4  IFIRV.LT.C. 130  TO  5 

DZ-D7/1L. 

GO  TO  3 
5 Z=Z*DZ 

GO  T0  3 

0 Z=V 

7 RETURN 

END 

3r0RtIS  ISOTAN 

SUBROUTINE  ISOTANI W.XMTt  XNT.CW.KWI 

FUNC(  W ) = W**XNT«-<  3 • ♦ < X H T ♦ 1 .)  / 7 ) * W * * ( X N T ♦ X MT  ) - C W 

L -D 

kw-f 

XCW-ALCG1L  ( 3 • * ( XMTf I . 1/7*1 
RM-W/(1C.**3. ) 

DW=W/1CD. 

0-W 

P3-FUNC (Q  1 

IF! A9SI RQI.LZ.RM 130  TO  1«4 
IF ( X NT  +XMT ritlC. 2 

1 1 = 1 

GO  T0  3 

2 I=-l 

3 IFCRO.GT.DIGC  TO  7 

<4  RW=FUNC(WI 

IF(A3S(RWI.LE.RMIG0  TO  1A 
KW=KW*1 

IF  (KW.GT.  100  ICO  TC  1*4 
IF  CRW.CT. 0 ICO  TC  5 
L = L + 1 

W=W-I*L*DW 

wv=w 

CALL  SIZECKCXNT  tXMTtXCWtWV.KX ) 

IF( KX. GT.O ) GO  TO  12 
GO  TO  A 

E V = W + I*  D W 

WV=V 

CALL  SIZECK CXNTt XNTf XCWtWVfKX) 

IF (KX.GT.0  > G C TO  12 
RV  = F 'JNC  ( V ) 

IF ( ABS (RV  ) •LE.RMI GC  TO  11 
KW=KW+1 

IF ( KW  »GT • 1 DO  I GC  TO  1*4 
Ir(  RV.GT.D  GO  TO  6 
OWrDW/lC. 

GC  TO  5 

0 W = W ♦ I * D W 

wv=w 

CALL  SIZECKfXNT tXNTf XCW. WVtKX ) 

IFIKX.GT.O )GC  TO  12 
30  TC  5 

C "C’JTINC  WHERE  FUNGI W)  IS  PLUS  FOR  W = I 
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7 RW=FUNCCWI 

IFC  A9RC  RWI .LE.RM  >30  TO  1*1 
IF  CKW.GT.1001G0  TO  1*1 
KW  = KW«-1 

IF  CRW.LT.U ISO  TC  3 
L-  L*  1 

W=W-*I  *LOW 
WV  = W 

CALL  SI7ECKCXNT  tXMT.XCWt  UV .KX I 
IMKX.GT.CIGO  TC  12 
GO  70  7 

3 V=W-I*OW 

WV  = V 

CALL  SIZECK C XNTt XHT tXCW tWV« KX I 
IF  CKX.ST.O  ICO  TC  12 
RV-FUNC  t V ) 

IF CABS (RV I.LE.RMIGO  TO  II 
K!J -KW«- 1 

IF  CKW.GT.  IGO  IGO  TO  1*4 
IFCRV.LT.PIGO  TO  9 
DW=DW/10. 

GC  TO  3 

9 W=W-I*OW 

WV  = W 

CALL  SIZECKC XNT  .XMTtXCUt WVtKX  I 
IFCKX.GT.CIGO  TO  12 
GO  TO  3 

C ROUTINE  FOR  CXMT«-XNT1=0 
1C  CHECK  = CW-3.*CXMT-»1.  1/7. 

T"C CHECK. LF.U IGO  TO  13 
y=CHECK**Cl./XNTI 
GO  TC  1*1 

11  W = V 

30  TO  1*C 

12  KV=20C 
30  TC  1*4 

13  KW=3C0 

1*4  RETURN 

END 

SFOR.IS  GIT  CK 

SU9R0UT INE  S IZ ECK C XNT . XMT .XCU»WV»KXI 
KX-C 

ST7ES=XNT*AL0G1CCWV) 

SIZEK=  CXMT+XNT I *AL0C1CCWVI 
STZEL=XCW*SIZCM 
IT  CABS (SIZESI.CT.2C.  IGO  TC  1 
IFC  A >3  S C S I Z Z M I .GT  .20.  ICG  TO  1 
IF ( ADS CSIZELI.GT. 2H. IGO  TO  1 
GO  TO  7 

1 KX=1 

2 RETURN 
END 
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